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Abstract 

We have generalized the approach in of Dunajski, Mason and Tod ^U] an d 
established a 1-1 correspondence between a solution of the universal Whitham 
hierarchy |2.'i| and a twistor space. The twistor space consists of a complex 
surface and a family of complex curves together with a meromorphic 2-form. 
The solution of the Whitham hierarchy is given by deforming the curve in the 
surface. By treating the family of algebraic curves in CP 1 x CP 1 as a twistor 
space, we were able to express the deformations of the isomonodromic spectral 
curve in terms of the deformations generated by the Whitham hierarchy. 

1 Introduction 

The Whitham hierarchy originates from the study of the 'dispersionless limit' 
of integrable systems. The typical setup involves an introduction of a small 
parameter e and a suitable rescaling of the time variables. For example, if one 
rescale the times x, t to x = ex, t — et in the KdV equation 

dtu — 6ud x u — d x u 

one obtains the KdV equation with small dispersion 

dtu(x, t, e) = 6ud x u(x, t, e) — e 2 d x u(x, t, e) (1) 

where d x and dt denote the derivatives with respect to the tilded variables. 

The problem is to study the behavior of the solution u(x, t, e) as e — > 0. In 
general, the solution contains parts that oscillates rapidly as e — > and the limit 
does not exist. 

However, the weak limit of u(x, t,e) can be studied. This is like 'averaging 
out' the oscillatory part and study the limit, and such an averaging process is 
called the Whitham averaging. 0, 0, [12, dU, [E], [H], E2, E3, 125], |23, 

M 

There is an interesting connection between the small dispersion KdV and 
deformations of Riemann surfaces. [T2] 
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A finite gap solution of the KdV equation can be expressed in terms of theta 
functions of a genus g algebraic curve (the spectral curve) . The spectral curve 
is fixed under the evolution of the finite gap solution. [TJ]], $ZU\ In many 
applications, a solution u(x,t) of the KdV equation can be approximated by a 
finite gap solution locally. However, the approximation becomes bad as (x, t) 
grows large and the spectral curve that classifies the finite gap solution varies 
slowly with respect to (x, t). To study the variations of the spectral curve, one 
could think of the spectral curve as depending on the 'slow times' (x = ex, t = et) 
and study the dependence of u on (x, t) in as e — > 0. 

In this case, The weak limit of u{x, t, e) can be described by 2g + 1 functions 
that satisfy the g-phase Whitham equations |2U 

d t Ui-Vg t i(ui,...,U2g+i)d x Ui = 0, 1 = 1, 2, . . . , 2g + 1 (2) 

where v 3i i{u\, . . . , u^g+x) depends on Uj on complete hyperelliptic integrals of 
genus g. These equations can in fact be written in a more compact form 

d t dp = d x dq (3) 

where dp and dq are certain meromorphic 1-forms that are normalized on the 
spectral curve with respect to a choice of a canonical basis of cycles. To be 
precise, the spectral curve in this case is a hyperelliptic curve 
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2g+l 

r(z) = n (* 



and dp, dq are 1-forms of the form 

dp 

dq 
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dp = / dq = 

>"2i J r 2i 

To obtain (j2J) from @, we expand dp and dq in terms of a z and coefficients 
would satisfy ©. [U, [12, 0, EH 

Krichever has studied equations of the type © on more general algebraic 
curves that cover the dispersionless limit of other integrable systems and in- 
troduced the notion of the universal Whitham hierarchy ■ The universal 
Whitham hierarchies, like are expressed in terms of meromorphic and holo- 
morphic 1-forms on Riemann surfaces. It has many interesting applications to 
2-D topological quantum field theory, Frobenius manifolds and string theory. 

n, 0, 0, m, eg, pa, esi 

In the case where the curve in the question is of genus 0, Dunajski, Mason 
and Tod has studied the dispersionless KP (dKP) equation from a point of view 
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of twistor theory, in which a solutions of the dKP equation is described by a 
family of rational curves in a complex surface and a meromorphic 2-form on 
the surface. This reveals interesting relations between the dKP equation and 
the Einstein- Weyl metric. [5], [3], |10|- In special solutions of the dKP 
equation was constructed by using twistor methods. 

In this paper, we have generalized the construction of JU| to the case where 
the curves are of arbitrary genus. We have established a one-one correspondence 
between a solution of the universal Whitham hierarchy and a twistor space, 
which is a complex surface with a family of genus g curves in it. To be precise, 
we have the following 

Definition 1 A twistor space T of the truncated universal Whitham hierarchy 
consists of: 

1. A 2- dimensional complex manifold T with a meromorphic 2-form II which 
is singular on a divisor D, 

2. A family of genus-g embedded curves {£ s ,t} and canonical basis of cycles 
on each curve, 

3. A covering Vp of a neighborhood U of the singular divisor D, and local 
coordinates fcj 1 on Vp such that k^ 1 has zeroes of order 1 at D. A fix point Pi 
on each curve S t such that Pi € D. 



Theorem 1 There is a one-one correspondence between a solution of the Whitham 
hierarchy and a twistor space T defined above. 

The idea is that as the curve moves around in the twistor space, the values of 
the 2-form II evaluated on different curves would give the values of the solution 
at different times. 

We have also studied the following problem which relates the isomonodromic 
spectral curve to the Whitham hierarchy. Consider the isomonodromic defor- 
mations of the system of linear ODE 



where A is a n x n matrix such that 



where G, F are n x r matrices constant in z, A = diag(ai, . . . , a r ) and C — 
diag(ci, . . . ,c„). 

Define the dual isomonodromic spectral curve |I5| . |28j of this isomonodromic 
system by 



These spectral curves are known to vary as the matrix A(z) deforms isomon- 
odromically. 



and 



A(z) = G T (A-z)- 1 F + C 
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The family of spectral curves is a family of algebraic curves in CP 1 x CP 1 
defined by polynomials P(z, w) that has degree at most n in w and r in z. Since 
the family of curves defined by all polynomials P(z, w) with degree at most n in 
w and r in z forms a twistor space for the Whitham hierarchy, we can compare 
the changes of the spectral curve induced by isomonodromic deformations and 
the changes induced by the Whitham hierarchy. This problem was posted by 
Takesaki in |3pJ, [3U| in a slightly different form. We were able to derive for- 
mulas that express the isomonodromic flows in terms of the Whitham flows in 
proposition |H] 

This paper is organized as follows. In section [21 we will give a brief intro- 
duction to the universal Whitham hierarchy and established some notations. In 
section [3] we will construct the twistor space for the Whitham hierarchy and 
in section 0] we will give some examples and study the isomonodromic spectral 
curves. 

2 The Universal Whitham hierarchy 

In this section we will give a brief introduction to the universal Whitham hier- 
archy in \2'6\ . 

Let M g ^M,n a be the moduli space of genus g curves Tg^ with N punctures 
P a (q), fixed n^-jets of local coordinates fc" 1 (q) in neighborhoods of P a (q) and 
canonical basis of cycles a, (q), bi(q). The marked points, local coordinates 
and basis change from curve to curve, but the genus g of the curves, the total 
number of marked points N and the numbers n a do not vary between curves. Let 
J2 a ( n « + 1) = H then the dimension of the moduli space M 9i N,n a is3g — 2 + H. 

That is, we have 

M g , N = {T g , v P a (q),k-\q) iai (q)M<l) e ^(Tg^Z)} 

The tautological bundle M g .N.n a —* Mg,N,n a is the fibre bundle over M Si jv inQ , 
with fibre at each point q of M g ^ tTla being the curve r 5i? . The curves T 9tq at 
different points are not holomorphic to each other in general. 

We consider each fibre r 9i9 as a genus g Riemann surface with canonical basis 
{cii (q) , bi (q)} of cycles, marked points P a (q) and local coordinates fc" 1 (q) near 
Pa (<?)■ We now define various meromorphic 1-forms dVL^ ^ (q) on these Riemann 
surfaces according to the following 

Definition 2 The meromorphic 1-forms c&Ia where A — (/i, v) a double index 
are defined by 

1. dQ aj i (q) are meromorphic 1-forms of the second kind which are holomor- 
phic outside P a (q) and have the form 

dn a , t (q) = d(kl l (q) + 0(k~ 1 (q))) (4) 
near P a (q), where i = 1, . . . ,n a . 
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2. d£l a fl (q), 1 are meromorphic 1-forms of the third kind with residues 
1 and -1 at P\ (q) and P a (q) respectively, and they look like 

dn afi = dk a ( q )(k- 1 ( q ) + o(k- 1 ( q ))) 

- dk 1 (q)(k- 1 ( q )+0(k- 1 (q))) (5) 

near P a (q) and Pi (q) . 

3. dQh,k{q), (k = l,...,g and h fixed) are the normalized holomorphic 
1-forms 

d> dn h , k (q) = <5 lifc , i,k = l,...,g (6) 

J ai(q) 

4- The differentials in 1. and 2. are uniquely determined by the normalization 
conditions 

I dn A (q) = (7) 

Ja,i{q) 

where the index A — (a, i). 

We now fix an index (a,i) = (1, 1) and denote dQi t i (q) by dp(q). Let the 
multi- valued function p (q) be denned by the abelian integral of dp (q) 

p{P,q)= [ dp(P,q) (8) 



where P G ^g,q- We will now suppress the dependence of q and bear in mind 
that everything depends on the point q of the moduli space. We can expand 
the local coordinates k a in terms of the function p near each marked point P a . 

CO 

k a = ^2 Va,s{p-Pa) s , a 7^1 (9) 

s=-l 

where p a — p (P a ), and for a = 1 

oo 

ki =P + J2»hsP s (10) 

We can now use the parameters 

[p a = p(P a ),i'a,8,a=l,---,N,s = -l,...,n a } 
a s = p(E a ), dp(S s ) = 0, a = l,...,2 5 -2 (11) 

Uf = <t dp, i = l,...,g 



as coordinates on the moduli space M g; jv,n a • Note that we only need the coeffi- 
cients JVs up to s = n a as we are considering n Q -jets, which means that terms 
of order k l a , i > n a are all equivalent and does not make any difference. 
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Suppose we have some unknown functions i Mj „ on the moduli space M gi jv,n Q 
which are in 1-1 correspondence with the 1-forms dil^^. That is, the t^ iV are 
unknown functions of the coordinates Also, suppose these functions 

form a coordinate system on a subspace M' of M gt N,n a ■ Let f2 M)1/ be the abelian 
integrals of <Kl^ M as in (JSJ. Let M' be the tautological bundle over M' . We 
now construct a 2-form IT on M', by using the and f^,„ as follows 

n = ^«l A Ae^ (12) 

A 

where for simplicity, we write (/i, v) — A and 

(5fi a = d p Vt A dp + ds^AdtB 
B 

When performing the derivative, we treat i A and p as independent variables on 
M'. 

Krichever's construction of the Whitham hierarchy is then equivalent to the 
following. One wants to ask: how should the tA depend on 111(1 if we want II 
to be simple? 

IlAn = (13) 

The coefficients of dp A dtA A dt B A die of II A II give the following set of 
partial differential equations 

^^^ABcdA^Bdp^lc = 

where the summation is taken over all the permutations of A, B and C and 
£abc is the sign of the permutation. By taking C = (1, 1), and denote t\ A by 
x, we get the following 

d A n B -d B n A + {n A ,n B } = (14) 

where {/, g} is the Poisson bracket with respect to the symplectic form dx A dp. 

This set of partial differential equations is called the truncated universal 
Whitham hierarchy W(n a ). 

By expanding the both sides of these equations in terms of p, we obtain 
partial differential equations of the coordinates Ijllfl with respect to the t A - By 
solving these equations, we obtain the coordinates l|ll|l as functions of t A and 
also defines M' . Of course, the t A do not form a complete set of coordinate on 
the moduli space M 9 jv,n Q , so the truncated Whitham equations in fact give the 
coordinates of the moduli space in terms of the t A on the subspace M' where 
the t A becomes a complete set of coordinate. 

The equations H14fl should be understood as a system of PDE relating the 
t A and the coordinates Ijlll) , which enter the equations as the coefficients of the 
multi- valued 'functions' Q A expanding in terms of the multi- valued 'coordinate' 
p on the curve ^ g , q - Namely, we choose branches of the functions tt A and p such 
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that the expansions give the coordinates The equations (|14|l should only 

be understood locally. Having this interpretation in mind, we would now take 
a closer look of the consequence on <|T1|) . 

Proposition 2 The truncated universal Whitham hierarchy \14[ ) gives 

d A Uf = d x Uf, U, A = ^dn A , i=l,...,g (15) 

d A a s = d x n A (S s ) , s = 1, ... ,2g - 2 

after expanding the terms near a s and taking the b-periods, where the 1- forms 
are defined in definition^ 

Proof. Consider the 6-periods of the equations (|14|) . First note that since 
p i ► p ~\~ Uf after going around the cycle bi, the vector fields d A becomes d A i— > 
9 A + d A Ufd p after going around a fr^-cycle. The &i-period of (|14|1 is therefore 

(d A uf - d x uf) d p n B - (d B uf - d x uf) d p n A + d A uf - d B uf (ie) 

By multiplying both sides by dp, we can view these as equations of 1-forms. 
Since the 1-forms df2 Mi „ are linearly independent, the coefficient of dQ B must 
vanish, which gives the first set of equations in (|22() . 

The second set of equations is obtained by looking at the behavior of ifhfjl 

i 

near the branch points cr s . Near the branch points er s , we have to use (p — a s ) 2 
as a local coordinate. We expand the functions fl A near a s as follows 

n A = Q A (E s ) + P A>1 (T) (p - os) 1 * + At2 (T) (p - a.) + ■ ■ ■ (17) 
We now expand l|14|) near the point a s . The term d A Q B is of the form 

d A n B = -d A a s d p n B + o (i) (is) 

and the term d x VL A d p VL B is of the form 

d x Q. A d p Vt B = d x VL A (a s )d p VL B 

+ Pa,! (T)0b,i (T) (-^j d x a s (p-ctsY 1 

- (^pA,i{T)f3 B fl{T)d x a s (19) 

+ \p A ,2{T)f3 B ,i{T)d x a s ) (p-a s )-i 
+ 0(1) 

since the second and the third terms in the right hand side is symmetric in A 
and B, the expansion of 114|) near a s is then 

(-d A a s + d x n A (e s )) d p n B + (-d B a s + d x n B (£.)) d p n A + o (i) (20) 
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we now choose fis such that <KIb (S s ) 7^ 0. For example, the holomorphic 
1-forms dQh.k only has 2g — 2 zeros while the number of the points E s is 2g. 
Therefore there exists some B = (h,k) such that <KIb (S s ) 7^ 0. By similar 
argument as before, we arrive at the second set of equations. □ 

We also want the truncated universal Whitham hierarchy W(n a ) to be em- 
bedded in an infinite hierarchy with n a — > 00. This would give us the following 

Proposition 3 As n a — > 00 for all a, the truncated universal Whitham hier- 
archy gives 

d A k a -{k a ,Q A }^0 (21) 
for the 1-forms VL A defined in ffil. 

Proof. The proof can be found in [231 - Let B — ( a ij) an d A = (f3,i) in 1)14(1. 
Let X- be the holomorphic part of X at P a . Then we have 

dAki - {ki,n A } = d B n A - {n A , (n B )_} 

Since A = ((3,i) and (CIb)- is holomorphic at P a , the worst pole that the right 
hand side could have at P a is of order i. Therefore 

d A ki-{ki,n A } = o(k^) 

By letting j — > 00, the proposition is proved. □ 

We can now define the truncated universal Whitham hierarchy as a system 
of PDE for the coordinates l(TT|) 

Definition 3 (Truncated universal Whitham hierarchy) The truncated uni- 
versal Whitham hierarchy W{n a ) is defined as the following system of PDE 

d A k a + {n A ,k Q } = 

d A Uf = d x U t A , U? = j dil A , i=l,...,g (22) 
d A a s = d x n A (S s ) , s=l,...,2g-2 
for the 1-forms defined in definition^ 

To define these equations, we only use the coefficients of the ft A up to the p na 
term. Since these coefficients are expressible in terms of 1)11(1 . the PDE 1(22(1 
defines the functional relation between the t A and 1(11(1 . 

The set of differential equations l(14|) with indices A, B = (a, i) only (that 
is, excluding the indices [h, k) in 4.) are called the basic Whitham hierarchy, 
they represents the dispersionless limits of integrable hierarchies. 

For simplicity, we would, from now on, call the truncated Whitham hierarchy 
W(n a ) the Whitham hierarchy. 
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Remark Although the vector fields 8a are multi- valued on M g ,N,n a , they are, 
nevertheless, well-defined as normal vectors of the curve T g q C M 9t N,n a - This 
turns out to be crucial in our construction of the twistor space. 



2.1 Generating form of the Whitham hierarchy 

The equations l|14|l can be considered as compatibility conditions of the system 

d A E = {E ,n A } (23) 

where Ep is defined on some open set Vp on T 9tq . We can perform another change 
of variable and view Ep,tA as independent variables. The set of equations l|14fl 
then become 

d A n B (E p ,t) = d B n A (Ep,t) (24) 

in these new coordinates. The equations (|24() suggest the existence of a potential 
S (Ep, i) such that 

d A Sp{Ep,t) = Q A {Ep,t) (25) 

The 2-form II is then 

II = 5Ep A SQp (26) 

where Qp = d E(i Sp. 

The functions Qp and Ep satisfies a 'string equation' 

{Qp,Ep} = \ (27) 

this implies that Qp is also a solution of (|23() 

d A Q0 = {Qf3,n A } (28) 



3 The twistor space of the Whitham hierarchy 
( genus- g curves) 

We will now proceed to construct the twistor space of the Whitham hierarchy. 
We will first construct a twistor space from the Whitham hierarchy and study its 
properties. We will then extract these properties to define a twistor space and 
recover the Whitham hierarchy from it. In constructing the twistor space (the 
direct construction), we will assume that we are given a solution to the universal 
Whitham hierarchy and therefore we know the functional relation between tA 
and the coordinates fTTJ) and also the subspace M' of M 9: M,n a where t A forms 
a complete set of coordinates. We will denote the tautological bundle over M' 
by M 1 . In particular, the 2-form II is defined on M' . 
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3.1 The symplectic reduction of M' 

Following the spirit of JU|> we wn l treat the 2-form n as a presymplectic form 
on M' and take the symplectic reduction of (M',Il). This then give us a 
2-dimensional manifold which we will call the twistor space of the Whitham 
hierarchy. 

We first study the kernel of the 2-form IT. Due to the simplicity condition 
(|13|) . the kernel is generically of codimension 2. We have the following 

Proposition 4 The kernel of the 2-form II is spanned by the vector fields 

d A - d x il A d p + d p 9, A d x (29) 

Proof. This can be verified by direct calculation. We first contract d A with II, 
where we use (p,t A ) as our independent variables. 

d A JU— {d A £lB — 9b^ A ) dts — d p fl A dp — d x fl A dx 

We now apply the (|14|l to the above equation and get 

d A JU— \_. (9x^b9 p CI a — dxtijidpfls) dig — d p Cl A dp — d x ti A dx (30) 

We then consider the contraction between II and d p , d x respectively. 
dp J II = dx + dpVLsdtB 

B 

d x Jli = -dp + ^d x n B dt B (31) 

B 

By comparing these with 1|30|) , we see at once that 

(d A - d x il A dp + d p il A d x ) J II = 

This concludes the proof of the theorem. □ 

We can now consider the distribution Y spanned by the vector fields (|29|) . 
Since these span the kernel of a closed 2-form, this distribution is integrable. 
Note that although the 2-form II is degenerate and singular on some codimension- 
1 sets, this distribution is still well-defined on these sets. We would like to 
consider the leaf space of this distribution. That is, we would like to study the 
space M' /Y . This will be our twistor space T. 

Remark Note that the distribution Y is well-defined despite the fact that 
both d A and £l A are not single- valued. The vector field d A — d x Vt A d p + d p Vl A d x 
is defined up to an addition of 

d A u 2 A d p - d x ufd p - d x ufd p n A d p + d p n A d x u?dp 
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which vanishes due to the second set of equations in definition [3J The distribu- 
tion is therefore well-defined. 

We can now study the properties of T. The following gives us a more concrete 
picture of how T looks like. 

Proposition 5 The variables Ep, Qp andk a are functions onT . In particular, 
if we use them as local coordinates on T , the 2-form II then becomes II = 
dEp A dQp on T. 

Proof. The proof is nothing more than unwinding the definition of {/, g} in 
equations (j2HJl, 1128(1 and 1)22(1 . which says that given a vector field X G Y, 
X(E)=X (Q) = X(k a ) = 0. □ 

There is also a family of embedded curves in T, which are projections of the 
curves ^ g , q - These curves are labelled by the parameter tA- In explicit terms, 
these embeddings are given by 

Q/3 = d E ,Sp (32) 

as this equation is defined on the curves T g , q for tA = constant. 

The 2-form II descends to the twistor space T to give a meromorphic sec- 
tion II of the canonical bundle of T, which is non-degenerate except on a 
codimcnsion-1 subspace. Given a section T of the normal bundle N q of a em- 
bedded curve r g g C T, the map /it = T J II|r gives a meromorphic 1-form on 
Tg iq which can only have poles at the singular set of S of II with order less than 
or equal to that of IT. This then defines an isomorphism between H° (T g<q , N q ) 
and H (T g , q , K <£> [D]) where K is the canonical bundle of the T giq and [D] 
is the line bundle associated to the pole divisor of II. Therefore we have the 
following 

Proposition 6 The normal bundle of an embedded curve T g ^ q S T is isomor- 
phic to K ® [D], where K is the canonical bundle ofT g q and [D] is the pole 
divisor of II. 

In fact, the 1-forms 8a-I II are nothing but the 1-forms — dflA, since we have 
^Jllls = (d x n A d P - d p n A d x ) J (cLEAdQ)| s 
= | (d x Qd p Q,A ~ d p Qd x VL A ) d p E 
- (d x Ed p fl A - d p Ed x Vt A ) d p Q }d P 
as dtA\n — 0. We now apply the string equation 127|) to obtain 

3 4 -l n| s = d p n A dp = -dfl A (33) 

3.2 The twistor space of the Whitham hierarchy 

This concludes our study of the properties of T. We will now define a twistor 
space independently and show that it gives a solution to the Whitham hierarchy. 
We first define the twistor space of a Whitham hierarchy 
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Definition 4 A twistor space T of the truncated universal Whitham hierarchy 
consists of: 

1. A 2- dimensional complex manifold T with a meromorphic 2-form II which 
is singular on a divisor D, 

2. A family of genus- g embedded curves {£ s ,t} and canonical basis of cycles 
on each curve, 

3. A covering Vp of a neighborhood U of the singular divisor D, and local 
coordinates fcj 1 on Vp such that k^ 1 has zeroes of order 1 at D. A fix point P\ 
on each curve S t such that P\ 6 D. 

The main result in this section is the following 

Theorem 7 There is a one-one correspondence between a solution of the Whitham 
hierarchy and a twistor space T defined above. 

Proof. Given a Whitham hierarchy, we can take the symplectic reduction of its 
moduli space as we did in the last section to get a twistor space T. The local 
coordinates k a near each marked point P a is well-defined on the twistor space 
due to f2"T|l and they give the local coordinates in 3. of definition 0] 

This gives the first half of the proof. To go the other way round, we want 
to recover the space M' and show that the pull-back II of the 2-form II to M' 
has the correct form. The simplicity condition follows automatically as II is the 
pull-back of a 2-form on a 2-dimensional manifold. 

The space M' is the set of curves {S ffjt }. By considering E and Q as functions 
of the tA, we can pull back the 2-form II to the space M g ^. This can be achieved 
by expanding E and Q near the the divisor D in terms of the local coordinate 




and consider the coefficients as functions of t. 

Now choose representatives 8a of normal vectors such that 

8 a ju = (d x n A d p - d p n A d x ) jn (34) 

where d x J II = —dp. We see that 

(d x £l A d p - d p ft A d x ) J II| S = d p Vt A dp 

since (Xg^ — —dp. We see that tA is a complete set of coordinates on M 1 because 
the 8a generates all the possible deformations of the curves and are independent. 
We have 

& = fe d A Edt A + dpEdp^j A d B Qdt B + dpQdp^j (35) 

We now contract II with Oa and make use of l)34|l . The dts component of 8a- 1 A 
then gives 

d x fldpEd B Q - d p Qd x n A d B E + d p n A d x Ed B Q - d p n A d x Qd B E (36) 
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while the dp component gives 

d p Q A d x Ed p Q - d p n A d x Qd p E (37) 

since this is just —dpiljy, we obtain the string equation 1|27|) . Note that l|34|) can 
be written as 

(d A - d x tt A d p + d p fl A d x ) (E) dQ = (d A - d x ft A d p + d p ft A d x ) (Q) dE 

this gives <|23[) and l|28[l . which implies the compatibility of the flows (|14f) . By 
applying ill'-SI) and i|28|) to (|36|l and making use of the string equation (|27|l . the 
expression <|3(il) becomes 

d x Sl A d p Sl B - d p VL A d x VL B = d B ft A - d A il B (38) 

this proves theorem follows. □ 

Remark As in the case of dKP equation JUj , the Whitham hierarchy depends 
on the choice of coordinates, which is reflected in 3. of definition Q] 

4 Examples 

We will now look at a few examples of twistor spaces of Whitham hierarchies. 
Example 1 The dKP equation] 1 (If 

The first example is the dKP equation, its twistor space was first constructed 
in JHj from the Einstein- Weyl metric. This is the case of a Whitham hierarchy 
with genus-0 curves. 

The twistor space of this equation is given by a 2-dimensional complex man- 
ifold T such that the 2-form II is singular at a connected set D with order 4. 
The Darboux coordinates E, Q are chosen so that E, Q are singular on D with 
orders 1 and 2 respectively. There is only one fixed point E = Q = oo. The 
local coordinate is chosen to be Q. The normal bundles of the curves are then 
O (2). (since the order of poles of the normal bundle at p — oo cannot exceed 
the correspondent one of Q). 

The solution to the Whitham hierarchy is recovered by expanding E and Q 
near z = oo. 

Q = p + J2 U *P~ 1 » E = ^2. Vl Q- % + x + Qy + Q 2 t (39) 

By using II = dEAdQ and the singular structure of II, we see that the pull-back 
II has the form 

fl = dx A dp + dy A d Qp 2 + u^j + dt A d Qp 3 + pu x + wA (40) 

the simplicity of this 2-form then gives a solution to the dKP equation 

{( Ul ) t - HI (ui) x ) x = ( Ul ) yy (41) 
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Example 2 Elliptic curves 

We now look at the set of elliptic curves T q in CP 1 x CP 1 , these are given the 
equation 

w 2 =4z 3 - gi z-g 2 (42) 

We take the twistor space to be CP 1 x CP 1 with coordinates (w, z) and consider 
the deformations of the curve. 

We first choose the 2-form to be dw A dz, and we will choose the local 
coordinate k to be — . The marked point will be the point where dw A dz be 
comes singular on the curve, that is, the point (oo, oo). We shall simply denote 
this point by oo. Let A be the uniformization parameter of the curve, and p (A) 
be the Weierstrass function. We choose our basis of cycle such that 

i_ i_ 

£d\ = 2g3 10 ^ 15 (43) 

Since z = p(A) and w — p' (A), we see that the divisor D is of order 7 on 
the curve. The space H° (T q , N q ) is therefore of dimension 7 (since there is no 
meromorphic 1-form with a simple pole on T q , the dimension is 7 instead of 
8). We also want to fix the marked point (oo, oo) so we want all the curves to 
intersect this point. This means we only consider a 6-dimensional subspace of 
H° (T q , Nq). In particular, since the normal vector N such that fipj has a pole 
of order 7 does not vanish at (oo, oo), it will move the marked point and we 
shall not consider flows generated by such normal vectors. 

We shall see in a moment that the full family of curves obtained from de- 
forming (|42|l is given by 

g 3 w 2 + g 4 w + g 5 wz+ = g 6 z 3 + g 7 z 2 - g x z - g 2 (44) 

subject to a scale invariant such that 

£M»=° (45) 

Each curve in (|44|) can be brought to the standard form (|42|l by a linear 
transformation 

w = aw + bz + c, z = ez + f 
(w'f = 4(z') 3 -~ 9l z'-g 2 (46) 

therefore we can write the Weierstrass functions as 

p(X)=ez + f, p'(A) = aw + bz + c (47) 

By comparing l|46|) with (|44|l . we see that the coefficients a, . . . , / are given by 

i 1 _i l_i 

a = g£, 0= -g 5 g 3 2 , c = -g 4 g 3 2 , 
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V4/ ' J lie 



.98 + -^gfas 1 



(48) 



91 



■7 ( .91 + 12e/ 2 - ^343553 1 



.92 = .92 - gif + if - c 

We can therefore compute various powers of the jet k in terms of the Weierstrass 
function and hence compute the 1-forms dflooj. We shall now denote df2oo,i by 



dn 2 = 
dn 3 = 

dflA = 



2e 



(p(X)-u)d\ 



4 GD (p'w--(pw-")) dx 

(£) 3 4p"(A) - ^p'(A) + [24/ + 6 (^) 2 ) (p(A) - W ) 



dA 



(;)" 



32p'(A)p(A) - — p"(A) 
e 



(49) 



24 



b\ 2 



+ 64/ p'(A) 



bf 

- 96— - 32c 

e 



6\3 



(p(A)-u) 



16(p'(A)) 2 + 16p(A)p"(A) - 160-p'(A)p(A) 

e 



+ (40(^) Z + 80/)p"(A) + 



A// 



320^- - 80c - 40 (^) 3 J p'( A ) 



^240(^) 2 / + 160^ - I651 + 10(^) 4 + 480/ 2 ^J (p(A) - w) 



(50) 



where u> is the the a-period of p(A)dA divided by the factor 
The holomorphic 1-form dClh 1, which we will denote by d£lo is 



dflr, 



2ff 3 10 3 7 15 



d\ 



63). 



(51) 



The Whitham hierarchy that we are trying to solve is the form 

djdSli - didSlj + {dCli, tf%} = (52) 

The coefficients of the expansion of the above equations can be expressed in 
terms of the <?, in (|44p. We will therefore treat the coefficients <7i as a solution 
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to l|52l) and express them in terms of U. In fact, we will express U in terms of 
the 

We first look at the tangent bundle of the curve T q . By differentiating (|44|) , 
we see that the tangent bundle is spanned by 

3 = {2g 3 w + 34 + g5z)d z - (3z 2 g 6 + 2zg 7 - g x - g 5 w)d w (53) 

and the normal bundle of the curve T q is spanned by 

d zQ ] 

wd w - zd z ,d w ,zd w ,d z , 7 - -,— — — -> (54) 

{2g 3 w + g± + g 5 z) (2g 3 w + g 4 + g 5 z) J 

since the contraction of these vector fields with dw A dz gives 1-forms of desired 
order. 

We can also see this as follows. Since the poles of dw, dz are of order 4 and 
3 at oo, d w and d z vanish to the order 4 and 3 at oo respectively. Therefore the 
vector fields wd w etc. are holomorphic on T q . However, the vector field w 2 d w 
and z 2 d z which are holomorphic on CP 1 x CP 1 are not holomorphic on T q . On 
the other hand, vector fields like zd w that are not holomorphic on CP 1 x CP 1 
becomes holomorphic on T q as T q only intersects z — oo and w = oo at (oo, oo). 

Also, the vector fields 775 p^—, r and 775 t^-. r are holomorphic as d w 

is tangent to T q when the denomination vanishes on T q , which we can see by 
setting 3 = in (|53|l . Therefore the above are holomorphic normal bundles on 

Since dm cicts on CP 1 x CP 1 by an infinitesmal translation of w, we see that 
the normal vectors H54J) generates the following movements of the curve. 

wd w - zd z = 5g 3 d g3 + 2g 4 d gi + 3g 2 d g2 

+ig 1 d gi + g 7 d 97 + 3g 5 d g5 
d w = -2g 3 d gi + g 4 d g2 + g 5 d gi , 

zd w = 2g 3 d 95 + g 4 d gi - g 5 d g7 , 

d z = g5d 9i + 3g 6 d g7 - 2g 7 d gi + gid g2 , 

= d n 



(2g 3 w+g 4 +g 5 z) u 92> 
da 



(2g 3 w+gi+g 5 z) "Si 

where we have used the scale invariance J2i 9$gi = to eliminate the ged ge 
term in the first equation. We see that the full family of curves is of the form 
l|44|) . By contracting the normal vectors with respect to the 2- form n and using 
Ha — —dQi, we see that 

d gi = — d tl +d gi t a d to 

96 



d 9 2 = gs'Qe'd, 



to 



d 93 — -p-d t5 + 7%5 t4 - 7^-9*3 + d g3 t d to 
5ff6 4 g 6 6g 6 
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2.9, 

.93 



■dt 2 + dg 4 todt 



(55) 



95 



.97 



*f6 J »6 z f6 

d 97 = --^9 t3 - TT'i^ti + "4^ti + 9g 5 tod to 
6 9e A 9e 9e 

while the functions <9 Si to are expressed in terms of gi and w. From 155(1 . we see 
that the functions 



h 
t 2 

U 
ts 



+ Cl 



2.9i + gf 

54 ~ 5557 
2^6 

5853 + gj 

5355 . 

45e 
5 3 2 



C2 



C3 



(56) 



556 



+ c 5 



where q are integration constants. While to satisfies the following equations 

-i 1 1 57 . 5 2 



d gi t 
dg 2 to 

dg 3 t 



dgjo 
dg 5 t 



_ 3 _ 2 

5 3 5 5 6 5 



3 \ .96 4g 6 g 3 



_j3 _ 2 

5 3 5 5 6 5 



1 

553 

5 2 5i 
2.9356 



,53 ' 2 6.96.93 3.96 4g$g 6 
K../- + ( — -4 91 



5455 

53 



5751 _ 2g| _ ff| _ gfg7 3g 5 g 4 
65653 3g 6 ^93 9e 5653 2g 3 

2~ ( _ 54<9 gi t ~ 9bd g . 2 t^j 

1 
253" 



5g54 555i 



65653 3g 6 



dg 2 t() 



( 9 ^ + 29597 - Ad t 

\ 65553 3.96 ' 4 y 91 ° 



d.97 f - 



1 

356 



5554 
253 



9ijd g2 to+(^ + 2g T \d gi t 



(57) 



17 



where e = (4 1 g 6 ) § . 

This gives a solution to the Whitham hierarchy (|52|l given by twistor theory. 

Example 3 Algebraic curves in CP 1 x CP 1 

This example generalized the last example to the case where the curves are given 
by a polynomial 

P(w,z) = 

In this case we choose the symplectic form to be dz A dw. The twistor space 
is then a subset of the normal bundle which consists of normal vectors T such 
that ht\t, are normalized 1-forms. This gives us the deformations of the curve 
£ that correspond to the Whitham deformations. 

Spectral curves of the dual isomonodromic deformations The isomon- 
odromic problem was studied by Jimbo, Miwa and Ueno. |16) . |17j Suppose we 
have an linear system of ODE 

^*(*) = A(z)V(z) (58) 
where A{z) is a n x n matrix-valued rational function in z 

where A a (z) are polynomials in z — a of degrees r a — 1, A oc (z) is a polynomial 
of degree and I? is a divisor in CP 1 . We call the system non- resonant if 
the leading terms of A a has distinct eigenvalues when r a > 1 and that no 
eigenvalues of the leading terms of A a differ by an integer when r a = 1. In this 
case, we can find unique formal solutions to the system of the form 

V(z-a) = r (Q) (z)exp(£> Q (z) + m a log(z — a)) (59) 
as z — > a, where Y^(z) = J2T=o ^ \ z ~ a )~ j is a 

formal Taylor series of ma- 
trix that is invertible at z = a and D a {z) is a diagonal matrix with polynomial 
entries in z — a. The matrix m a is a diagonal matrix and is called the exponent 
of formal monodromy. 

It is known (e.g. 122) that near each pole a of A(z), there exist sectors Sa 
on which there exist solutions that are asymptotic to the one in (|59|l . 

By comparing solutions on different sectors, one obtains the Stokes matrices 

Sjlp = y&cPj 't&aP ■ Similarly, analytic continuation of solutions near differ- 
ent poles define the connection matrices. JHI Let us denote the coefficients of 
D a and the pole positions a by T collectively. The isomonodromic problem is to 
find the dependence of the matrix A(z) on T such that the Stokes matrices, the 
connection matrices and the exponent of formal monodromy remain constant. 
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It was shown in that the dependence A(z,T) of A(z) on T is determined 
by the following differential equation of A(z, T) 



dA(z, T) = d z B(z, T) - [A{z, T),B(z, T)} (60) 

where d is the exterior derivative with respect to T and B{z 1 T) is a rational 
matrix 1-form. 

The 'spectral curve' of the isomonodromic problem is then defined by the 
algebraic equation EB| 

det (w — A{z, T)) = (61) 

Unlike the spectral curves of other integrable systems, the 'spectral curve' of an 
isomonodromic problem is not preserved by the isomonodromic flows. One can 
see this from Pty. The term [A(z,T),B(z,T)\ on the right hand side of JBDJl 
would preserve the spectral invariance whereas the term d z B(z,T) would lead 
to deformations of the spectral curve. 

From now on we will consider a simple case in which the matrix A(z) is of 
the following form 

-W-tfr^y + c < 62 > 

where Ai are rank 1 matrices constant in z and C is a diagonal matrix constant 
in z, C = diag(ci, . . . , c„). The deformation parameters in this case are the 
and the Cj. 

In this technique called 'duality' ^S] can be employed to obtain an 

algebraic curve that is birational to the spectral curve. Moreover, this curve is 
defined by the zero locus of a polynomial P(z, w) = in CP^CP 1 . These curves 
were considered by Sanguinetti and Woodhouse. [2E| The explicit construction 
goes as follows. Let A be the diagonal matrix A = diag(ai, . . . , a r ). Since all 
the Ai are of rank 1, we can write the matrix A(z) as follows 

A(z) = G T {A-z)- 1 F + C (63) 

where G, F are n x r matrices constant in z. 

There is a 'dual isomondromic problem' associated to it, which is defined by 

d w ~ F{c~~w)- 1 G T -A (64) 

It was shown |15| that when A(z) is deformed isomondromically, the dual system 
also undergoes isomonodromic deformation. 

The dual isomonodromic spectral curve is defined by the following polynomial 
in w and z 



det 



M - 




°)1 




(o 


W J 



= (65) 
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where M is the following matrix 



A —F 
G* G 



The genus of such curve is shown in [2Hj to be (n — l)(r — 1). The whole set 
of curves £t obtained by deforming the spectral curve through isomonodromic 
deformations form a n + r parameter family of curves. This is a subset of curves 
defined by the zero locus of all polynomials 

r n 
i=0 j=0 

which is a, nr — 2n + 2r + g parameter family of curves. We can treat this family 
of curves and the 2-dimensional complex manifold CP 1 x CP 1 as a twistor space 
of Whitham hierarchy as in example [3] The 2-form is given by dz A dw and 
the fixed points are the intersection points between z = oo, w = oo and the 
curves. These are the points (z — oo,w — c,-) and (z = cti,w — oo). The local 
coordinates near z = oo are chosen to be z~ x and the local coordinates near 
w = oo are chosen to be w . In this case there is no natural choice of the point 
Pi and we will choose it to be [z = cxi,w = oo). 

In this case the Whitham times can be solved by the hodograph method. 
|25) . EH They are given by 

tcti.j — ies w=00 w~ :l zdw, j = 0,1 

t Cj j = — res z=00 z~ 3 wdz, j — 0,1 (66) 

th,j = f zdw, j = l,...,g 



where we used the index on to denote the point (z = on, w — oo) and the 
index a to denote the point [z = oo,w — c,). The residues are taken around 
the corresponding points. The integrals defining the th t j are around a certain 
choice of a-cycles on the curve. Note that these formula are independent on the 
choice of Pi . 

Since the Whitham flows generates all the deformations of the algebraic 
curve, we can express the isomonodromic deformations of the spectral curve in 
terms of the Whitham flows. 

d = dtA g 

To see how the isomonodromic times are related to the Whitham times, we first 
consider the behavior of w near (z — oo,w = Cj). Near these points, w has the 
Laurent expansion 

w = c t + m^z- 1 + 0(z- 2 ) (68) 
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where is the i th entry in the formal exponent of monodromy at z = oo. 
From 166(1 we see that 

tcifi = -W'ooi *Ci,l = -c. (69) 

Similarly, near (z = a^,w = oo), z has the expansion 

2 = on + n^w' 1 + 0(w~ 2 ) (70) 

where is the i entry in the formal exponent of monodromy at w = oo of 
the dual system ijBlfr. We see that 

*a S :0 = ^lo) *on,l = a i (71) 

By substituting these into 1|67[) we obtain 

■S— = d ai ,i + V* -S— ( / ^rfw ) d h k 



d_ 



9c i< i+J2j-( ( f zdw ) d h,k (72) 

k— 1 ^ \J a k / 

To calculate the second terms, first note that since z is single value on the spec- 
tral curve, the derivatives (jf aM zdivj and [J> ak zdw\ has no boundary 



terms. That is, we have 

S- ( & zdw] = <fc -^-zdw (73) 
dci \J ak J J ak dCi 

( I zdw] = I zdw (74) 
dai \J ak J J ak dan 

To evaluate d Ci zdw or d ai zdw, we can either think of the isomonodromic flows 
as flows that keep z fixes and changes w as in the original isomonodromic system 
<|58[) , or we could think of these flows as keeping w fix and changing z as in the 
dual isomonodromic system (|64|) . To see that these two different interpretations 
make no difference in the final result, let us denote the determinant in l|65|l by 
P(z,w,T). Let 

r n 

P(z,w,T) = ^^X ij (T)z i wi =0 

Then we see that 

(r n \ 
^^^.i.^jzV {d,p{z,w,T)y x 
»=0 3=0 J 

(r n \ 
^^^.^(rjzv (a w p(z, w ,T))- 1 
i=0 j=0 / 
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and similar expressions for d Ci w and d Ci z. If we think of the isomonodromic 
flows as flows that fix w, the integrals in l|73(l are then 

d zdw = - & [ yy^I.jfTjz'^ I (dzPiz^w^T))- 1 dw (75) 



da 



on the other hand, if we think of the flows as flows that fix z, the integrals 
become 



-^-zdw=<j) ( y^y^a Ql X y (T)zV \{d w P(z,w,T))~ 1 dz (76) 

since £ zdw — — £ wdz as zdw + wdz is a total differential. 

k 

To see that the right hand sides of (|75|l and l|76|) are the same, note that 

by differentiating the expression P(z, w, T) = 0, we have (d w P(z, w, T)) dz — 

— (d z P(z,w,T))~ dw. By replacing the on derivatives in (|75|) and lj7T)l) by 

derivatives of Cj and apply similar argument, we see that the same holds for the 

Ci derivatives. 

To compute the derivatives of the Xij , recall that an isomonodromic defor- 
mation changes the matrix A(z) in the following way 

dA(z, T) = d z B(z, T) - [A(z, T),B(z, T)] 

The matrix 1-form B(z 1 T) in this case is 

r n 

B{z, T) = J2 {z, T)da t + B c (*, T)dc t 



where the B ai (z,T) and the B Ci (z,T) are as follows JS] 

(z - on) 

B c .(z,T) = zE i + yy EiA <** E i + E i^ E i (78) 

j^ti k—l J 

where Ei are the n x n diagonal matrices with 1 on the i th entry and zero 
elsewhere. 

The key observation here is that d z B ai = d^ .A(z, T) and d z B Ci = d^.A(z, T), 
where d e a . and d e c . denotes the explicit derivatives with respect to cti and Cj. 
Therefore we have 

diA(z, T) = d*A(z, T) - [A(z, T), B,(z, T)} 

where the index i will be used to denote cti or &. 
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Since the commutator term [A(z, T), Bi(z, T)] fixes the spectral invariance 
of the matrix A(z, T) and it is known that the coefficients Xij(T) of the spectral 
curve are spectral invariance, |15| we see that the derivatives dxXij in l|7ti|l or 
(|75|l are the same as their explicit derivatives with respect to the times T. That 
is 

drXij — d^Xij 

this observation simplifies the calculations significantly as the explicit depen- 
dence of the coefficients Xij on a, and Cj are polynomials, which can be seen 
from l|65[). 

In summary, the isomonodromic flows and the Whitham flows are related 

by 

Proposition 8 Let A(z) be a matrix of the form 

A(z) = G T (A-z)- 1 F + C 

where G, F are n x r matrices constant in z, A = diag(«i, . . . , a r ) and C = 
diag(ci, . . . ,c„). 

Define the spectral curve of A{z) by 



dot 

where M is the following matrix 



M - 




°)1 




(o 


w J 



= 



A —F 
G* C 



then, when the system of linear ODE 

dY 

?f = A(z,T)Y 
az 

deforms isomonodromically with deformation parameters cxi and Ci, the defor- 
mations of the spectral curve in CP 1 x CP 1 can be expressed in terms of the 
Whitham flows as follows 

g „ / r n \ 

d m = d ai)1 +J2f ^^^^(rjzV \p- l dzd Kk 

k=l '' a k \ i=0 3=0 J 
9 r r n \ 

9 Ci = d CtA + J2f ^^^(TlzV P ffi Ht (79) 

fc=i -' a fc yi=o j=o J 

where df denotes the explicit derivative and P w is the w derivative of P(z, w, T). 

Here the twistor space of the Whitham hierarchy is taken to be the set of 
algebraic curves given by the zero locus of polynomials P(z,w) in CP 1 x CP 1 , 
where P(z, w) has degree at most n in w and r in z. The marked points are the 
points where z or w is infinity. 
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Remark As one could see from (I79f) that, the contraction diJH of di with 
the 2-form dz A dw would result in a 1-form on the curve that is holomorphic 
everywhere apart from a second order pole at the point (z = a.i,w = oo) for 
d ai or (z = oo,w = Ci) for d Ci . In POli Takasaki rose the question of 

whether it is possible to choose canonical basis of cycles that could make these 
1-forms normalized. This is the same as choosing a basis of cycles so that the 
summation parts in l|79|l vanish. We do not know the answer but these formula 
would simplify the problem significantly since all the explicit derivatives can be 
computed. 

We would now give an example. 

Example 4 Painleve V 

The Painleve V equation is the following second order nonlinear ODE 



u 



1 1 \ u' (u-lf ( P\ ju 8u(u + l) 



= («') Ur + — r -T + ^r^ + ~ 



2u u — 1 / t t 2 \ u j t u — 1 



where a, (3, 7 and 8 are constants. It can be represented as the isomonodromic 
deformations of a system of linear ODE in which A{z) is of the following form 

EH! 

A(z) = - — - 



z 



-t 



We will follow the approach in ^S] to construct the spectral curve. We first set 
the constant 5 to be 5 = 2. In this case, the matrices A® and A\ will be of rank 
1 and a factorization in (jfiHJ) can be carried out. 

In |15| , the Painleve V was represented by the isomonodromic deformations 
of the system of linear ODEs 

My-l = ft *\ J- Xiyi . 

dz \ -t J \ Ji 



2(z-l) 



-X2IJ2 - M2 

X% X 2 V2 




in which fi-y, 112 and c = h{xiyi + X2IJ2) are constants. These constants are 
related to a, (3 and 7 as follows 

«=f, /3 = -f, 7 = 4c + 2 

x 2 

The solution y to the Painleve V equation is given by u = — j. 
The factorization (|6;-S[) is given by 

F = j_{xi vi r = ±fyi + ^ ~ Xl 
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with A = diag(0, 1). 

The spectral curve is then 

P(z, w) — z 2 w 2 — w 2 z + (ni)wz — (fix + H2)w — t 2 z 2 + (t 2 + txiyi)z 

This curve is of genus 1. The fixed points are (z = 0, w = oo), [z = 1, w = oo), 
(z = oo, w = t) and [z = oo, w = —t). We can now apply H79|) to express dt in 
terms of the Whitham flows 

d t = <9 U + I (-2tz 2 + (2t + x iyi )z - 2c) P^dzdh,! 

J a 

where a is an a-cycle chosen on the curve. 
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